Proceedings The Art of Modelling Stars in the 21st century 

Proceedings IAU Symposium No. 252, 2008 © 2008 International Astronomical Union 

Licai Deng K.L. Chan & C. Chiosi, eds. DOI: 00.0000/X000000000000000X 



Thermohaline mixing in low-mass giants 

M. Cantiello 1 and N. Langer 1 

1 Astronomical Institute Utrecht, Princetonplein 5, 3584 CC Utrecht, The Netherlands 
email: M.Cantiello@uu.nl, N.Langer@uu.nl 

Abstract. Thermohaline mixing has recently been proposed to occur in low mass red giants, 
with large consequences for the chemical yields of low mass stars. We investigate the role of 
thermohaline mixing during the evolution of stars between 1 Mq and 3 Mq , in comparison to 
other mixing processes acting in these stars. We confirm that thermohaline mixing has the 
potential to destroy most of the 3 He which is produced earlier on the main sequence during the 
red giant stage. In our models we find that this process is working only in stars with initial mass 
M < 1.5 Mq. Moreover, we report that thermohaline mixing can be present during core helium 
burning and beyond in stars which still have a 3 He reservoir. While rotational and magnetic 
mixing is negligible compared to the thermohaline mixing in the relevant layers, the interaction of 
thermohaline motions with differential rotation and magnetic fields may be essential to establish 
the time scale of thermohaline mixing in red giants. 
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1. Introduction 

Thermohaline mixing is not usually considered as an important mixing process in single 
stars, since the ashes of thermonuclear fusion consists of heavier nuclei than its fuel, and 
stars usually burn from the inside out. The condition for thermoh aline mixing, however. 



is that the mean molecular weight (/i) decreases inward. Recently ICharbonnel k. Zahn 



(|2007al CZ07) identified thermohaline mixing as an important mixing process which 



significantly modifies the surface compositi on of red giants after th e first dredge-up. The 
work by CZ07 was triggered by the paper of Eggleton et al. ( 20061 EDL06), who found a 



/i-inversion in their 1 Mq stellar evolution model, occurring after the so-called luminosity 
bump on the red giant branch (RGB), which is produced after the first dredge-up, when 
the hydrogen-burning shell reaches the chemically homogeneous part of the e nvelope. 



Theu-inversion is produced by the reaction 3 He( 3 He,2p) 4 He (as predicted by lUlrich 



(1972)). R does not show up earlier, since the magnitude of the ^-inversion is small and 



negligible compared to a stabilizing /i-stratification. 

The mixing process below the convective envelope in models of low-mass stars turns 
out to be essential for the prediction of the chemical yield of 3 Hc (EDL06), and to under- 
stand the surface abundances of red giants, in particular the 12 C/ 13 C ratio, and the 7 Li, 
carbon and nitrogen abundances (CZ07). This may also be important for other occur- 
rences of thermohaline mixing in stars, i.e., in single stars when a /i-inversion is produced 
by off-center ignition in semi-degenerate cores, or in st ars which accrete chem ically en- 



riched matter from a companion in a close binary (e.g. JStancliffe et al. 1120071) . Accreted 
metal-rich matter during the phases of planetary formation also leads to thermohaline 
mixing. The host stars of exoplanets present a metallicity excess compared to stars in 
which no planets have been detected. This metallicity excess can be reconciled with the 
overabu ndances expecte d in cases of accretion if thermohaline mixing is included in the 



picture (jVauclair 1 120041 ) 
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2. Method 



We use a 1-D hydrodynamic stellar evolution code (jYoon et al. 1120061 and references 
therein). Mixing is treated as a diffusive process , the contributions to the diffusion coef- 
ficient are convection, semiconvection, thermohaline mixing, rotationally induced mixing 
and magnetic diffusion. The code includes the effect of centrifugal force on the stellar 
structure, and the transport of angular momentu m is also treated as a diffusive process 



(|Endal fe Sofia I Il978t IPinsonneault et al. I fl989h . The condition for the 



occurrence 



of 



thermohaline mixing is 



^ V - V ad < 



(2-1) 



i.e. the instability operates in regions that are stable against convection (according to the 
Lcdoux criterion) and where an inversion in the mean molecular weight is present. Here 
<p = (d\np/dlnp,) P>T , S = ~{d\np/d\nT) P ^ V M - dlnp/dlnP, V ad - (31nT/dlnP) ad , 
and V = dlnT/dlnP. Numerically, we trea t thermohaline mixing through a diffusion 
scheme (|Braun |[l997 ; Wellstein et al. I l200lh - Th e cor responding diffusion coeffi cient is 
based on the work of Stern (1960). IUlrich I (jl972h . and IKippenhahn et al. 1 (|l980h : 



i£\7 



3K 



2 pc P (V ad - V) 



(2.2) 



where p is the density, K — AacT 3 / (3np) the thermal conductivity, and cp = (dq/dT)p 
the specific heat capacity. The quantity a t h is an efficiency parameter for the thermoha- 
line mixing. The value of this parameter depends on t he geometry o f the fingers arisin. 



from the instability and is still a matter of debate (|Ulrich 1 11972c IKippenhahn et al 



1980t Charbonnel fc Zahn 2007a ) . Most of the calculations have been p erformed with 



2.0, corresponding to the prescription of IKippenhahn et al. (1980), although we 



also investigated the effect of using different values of a t h- 

In the code rotational mixing is included. Four different diffusion coefficients are cal- 
culated for dynamical shear, secular shear, Eddington-Sweet circulation and Goldreich- 
Schubert-Fricke instability. Details on the physics o f these instabilities and their imple- 



mentation in the code can be found in lHeger et al. I (|2000h 



Ch emical mixing and transport of angula r momentum due to magnetic fields ([Spruit 
2002) is included as in iHeger et al. 1 {2000). The contribution of magnetic fields to the 



mixing is also calculated and added to the total diffusion coefficient D entering the 
diffusion equation. 

We compute evolutionary models of 1.0 Mq, 1.5 M Q , 2.0 M© and 3.0 M at solar metal- 
licity (Z=0.02). The init ial equatorial v elocities of these models were chosen to be 10, 



45, 140 and 250 kms 1 ( Tassoul 2000l) : we assume the stars are rigidly rotating at the 



zero-age main s equence. Throughout the evolution of all models, the mass-loss rate of 
Reimers ( 19751) was used. 



3. Results 

We confirm the presence of an inversion in the mean molecular weight in the outer wing 
of the H-burning shell, after the luminosity bump on the red giant branch. According to 
inequality (|2.1j) this inversion gives rise to thermohaline mixing in the radiative buffer 
layer, the radiative region between the H-burning shell and the convective envelope. 
Results of our calculations for the 1.0 and 2.0 M Q models are shown in Fig.[TJ In our 1 M 
model, thermohaline mixing develops at the luminosity bump and transports chemical 
species in the radiative layer between the H-burning shell and the convective envelope. 
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Figure 1. Left: evolution of the internal structure of a 1.0 Mq star from the onset of ther- 
mohaline mixing to the AGB phase. Green hatched regions indicate convection, yellow filled 
regions represent semiconvection and red cross hatched regions indicate thermohaline mixing, 
as displayed in the legend. Blue shading shows regions of nuclear energy generation. Right: 
same as left panel, but for a 2.0 Mq star. 



This results in a change to the stellar surface abundances. The left panel of Fig. [2] shows 
the evolution of the 3 He surface abundance, and of the ratio 1213 at the surface as a 
function of time in the 1 M Q model, confirming the result of EDL06 and CZ07, namely 
that thermohaline mixing is efficiently depleting 3 He and lowering the ratio 1213 on the 
giant branch. Interestingly the 2M model shows a different behaviour: thermohaline 
mixing is not able to connect the H-burning shell and the convective envelope, resulting 
in no change to the surface abundance of 3 He and to the ratio 1213 after the luminosity 
bump. Our models with ath = 2.0 show thermohaline mixing to be important during the 
RGB only in stars with mass M <J 1.5 Mq. 

While CZ07 and EDL07 investigate thermohaline mixing only during the RGB phase, 
we followed the evolution of our models until the thermally pulsating AGB phase (TP- 
AGB). In fact a /i-inversion is always created if a H-burning shell is active in a chemically 
homogeneous layer; this happens not only during the RGB phase, but also during the 
horizontal branch (HB) and the AGB. The size of the /^-inversion depends on the local 
amount of 3 He. As a consequence, the efficiency of thermohaline mixing in different 
evolutionary phases depends on the amount of 3 He left by previous mixing episodes. 

After the core He-flash, helium is burned in the core, while a H-burning shell is active 
below the convective envelope. We found that during this phase thermohaline mixing is 
present and can spread through the whole radiative buffer layer in our 1 Mq model (left 
panel in Fig. [I). In this model the surface abundances change also during this phase 
because the H-burning shell and the envelope are connected. This is shown in Fig. [2l 
left panel, where surface abundances change after the l uminosity peak cor r espon ding 
to the He-flash. We stress that using the prescription of Kippenhahn et al. ( 1980h for 



thermohaline mixing allows our model to reach this phase without completely burning 
the 3 He; models of CZ07 almost completely depl ete 3 He in the envelope during the RGB 
phase because of their higher diffusion coefficient ( Ulrich 119721 ). In this case thermohaline 
mixing would be much less efficient, during the subsequent evolutionary phases, due to 
the lower abundance of 3 He. 

The last nuclear-burning phase of a low-mass star is the AGB, and is characterized 
by the presence of two burning shells and a degenerate core. The star burns H in a shell 
and the ashes of this process feed an underlying He shell. During the most luminous 
part of the AGB the He shell periodically experiences thermal pulses (TPs); in stars 
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Figure 2. Left: evolution of the surface abundance of the 12 C/ 13 C ratio (dotted red line) and 
3 He (dashed green line), and of the luminosity (solid blue line) from the onset of thermohaline 
mixing up to the AGB phase for a 1.0M Q star. Right: diffusion coefficients in the region 
between the H burning shell and the convective envelope for the 1.0 M© model during core 
He-burning (t = 1.267 x 10 10 years). The black, continuous line shows convective and thermo- 
haline mixing diffusion coefficients, the green, dashed line is the sum of the diffusion coefficients 
due to rotational instabilities while the blue, dot-dashed line shows the magnitude of magnetic 
diffusion coefficient. 



more massive than about 1.2 M© these TPs are associated with a deep penetration of 



the c onvective envelope, the so-called third dredge-up (3DUP) (jWallerstein fc Knapp 



1998). We find thermohaline mixing to be present during the TP- AGB phase. Depending 
on the mass of the model the diffusion process is able to connect the H-burning shell with 
the convective envelope during the whole interpulse phase. In a 1 M Q model thermoha- 
line mixing connects the H-burning shell to the convective envelope (Fig. [3J left panel) , 
confirming that this mixing process is more efficient at lower masses. The occurrence of 
thermohaline mixing in this late evolutionary phase is critically determined by the mix- 
ing history of the star. This is because the amount of 3 He present during the TP- AGB 
phase is controlled by the occurrence of thermohaline mixing and by its efficiency (ath) 
in previous evolutionary phases. 

Thermohaline mixing during the TP- AGB is a circulation of the type inferred by 
Cameron & Fowler (1971), since it carries 7 Be out of the high-temperature zone. 7 Be 
can then produce 7 Li by electron capture, which is transported in the convective enve- 
lope. This re sults in surface enrichm ent of 7 Li, as shown in the right panel of Fig. [3l 
Interestingly lUttenthaler et al. ( 2007 ) recently reported the detection of low-mass, Li- 
rich AGB stars in the galactic bulge. Given their low mass, these stars are not expected 
to experience any hot bottom burning. We argue that thermohaline mixing is a possible 
explanation for their abundance anomalies. 



4. Rotation and magnetic fields 

In our models we found that in the relevant layers thermohaline mixing has generally 
higher diffusion coefficients than rotational instabilities and magnetic diffusion. This 
result is not valid for magnetic stars, stars that p ossess anomalous surface fie lds of a few 
10 2 to about 10 4 G, believed to be of fossil origin (jCharbonnel fc Zahn Il2007bl ) . The right 
panel of Fig.[2]clearly shows that rotational and magnetic mixing are negligible compared 
to the thermohaline mixing in our 1.0 M© model. The only rotational instability acting 
on a shorter timescale is the dynamical shear instability, visible in the right panel of Fig. 
[2] (dashed line) as a spike present at the lower boundary of the convective envelope. This 
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Figure 3. Left: evolution of the region between the H burning shell source and the convective 
envelope during a thermal pulse in a 1.0 M© star. Green hatched regions indicate convection and 
red crossed regions indicate thermohaline mixing. Blue shading shows regions of nuclear energy 
generation. Right: evolution of Lithium surface abundance during one thermal pulse in a 3Mq 
model. The black, continuous line shows a model evolved with a t h = 1000; the blue, dotted line 
refers to the same model evolved with a thermohaline mixing efficiency ath =200. In both cases 
the model is experiencing third dredge-up. The evolution of the star prior to the TP-AGB has 
been calculated with ath — 2. 

instability works on the dynamical timescale in regions of a star characterized by a high 
degree of differential rotation. However, if present, this instability acts only in a very 
small region (in mass coordinate) at the bottom of the convective envelope. As a result 
thermohaline mixing is still setting the timescale for the diffusion of chemical species 
from the convective envelope to the hydrogen-burning shell. 

The interaction between rotation and thermohaline mixing is more difficult to address, 
since it requires full hydrodynamic calculations. We expect the speed of thermohaline 
mixing to be affected by differential rotation and rotational instabilities, since these are 
able to change the geometry of the fingers, and to create turbulence, respectively. Canuto 
et al. (2008) recently suggested that turbulence must be taken into account to explain 
both laboratory and ocean data of double diffusive processes. 

The interaction of magnetic fields with thermohaline mixing must also be considered. 
This depends on the geometry and the magnitude of the magnetic field: if the magnetic 
field is stronger or of the same order of the equipartition value the magnetic field al- 
lows plasma motions only along the direction of the field lines. This could result in the 
inhibition of thermohaline mixing, as discussed with detailed calculations by Zahn and 
Charbonnel (2008). However strong magnetic fields are observed in only a small fraction 
of low-mass stars (Wolff (1968); Power et al. (2007)). Concerning the interaction of weak 
fields with thermohaline motions, a change in the speed of the fingers is also expected 
and fully MHD calculations are needed to understand the process. 

5. Discussion 

We comfirm the results of EDL06 and CL07: thermohaline mixing in low-mass giants 
is capable of destroying large quantities of 3 He, as well as decreasing the ratio 1213. 
Thermohaline mixing starts when the hydrogen burning shell moves into the chemically 
homogeneous layers established by the first dredge-up. Our models show further that 
thermohaline mixing remains important during core helium burning, and can still be 
relevant during the AGB phase — including the termally-pulsing AGB stage. This can 
result in important changes to the surface abundances of low-mass stars. The quantitative 
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discussion is complicated by the fact that the capability of thermohalinc mixing to change 
surface abundances depends on an uncertain efficiency parameter, as well as on the local 
3 He abundance. The efficiency parameter a t h is still a matter of debate and is probably 
strongly affected by the interaction of thermohalinc mixing with rotation and magnetic 
fields. The local 3 He abundance depends on the previous history of mixing. 

Our calculations show that in the relevant layers, thermohalinc mixing generally has a 
higher diffusion coefficient than rotational instabilities and magnetic diffusion. Still, the 
interaction of thermohaline mixing with magneto-rotational instabilities is important: 
we expect the speed of the mixing to be strongly affected by the presence of differential 
rotation and magnetic fields. To better understand the picture it would be desirable to 
have realistic MHD simulations of thermohalinc mixing. 
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